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1. Introduction 



Exactly solvable lattice models have numerous applications to statistical physics, 
quantum field theory and solid state physics. Effective methods for calculating 
partition functions based on analytical properties were proposed by Baxter. 1 Re- 
cently, the Kyoto group ( Jimbo, Miwa, Nakayashiki, . . . ) considerably simplified 
these methods and developed an approach to calculating correlation functions at 
finite distances. 2-8 They applied this approach to the antiferromagnetic phases of 
the six-vertex 2 ' 3 and eight-vertex 7 ' 8 models and to the SOS and RSOS models. 4-6 
An interesting challenge is to generalize this approach to ferromagnetic and dis- 
ordered phases of the eight-vertex (and six-vertex) model. The problem is that 
the usual arguments for applicability of the corner transfer matrix approach break 
down in these phases. On the other hand, Baxter found some symmetries between 
different phases of the eight- vertex model. 1 More precisely, the partition function 
is invariant under some transformations of local wieghts, and these transforma- 
tions connect values of weights corresponding to different stable phases. We show 
that these symmetries can be extended to correlation functions. Then one can 
express correlation functions in ferromagnetic and disordered phases in terms of 
correlators in the antiferromagnetic phase. We formulate the rules for calculating 
correlation functions directly in the ferromagnetic and disordered phases using 
appropriate parametrizations of the local weights. 

Now we state the notations. 'Spins' of the eight-vertex model e = ± = ±1 
are situated at the links. The interaction is attached to the vertices and described 
by the weight matrix (or .R-matrix) R%\%\ (Fig. la) 
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The partition function is given by (see Fig. 2) 

z= e 



We need one more element (Fig. lb). It is the spin flop operator a 1 = ^ ^ ^ ^ . 
It may be considered as a kind of a two-link vertex. 
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Fig. 2 



There are five regions of parameters 1 a, 6, c, <i 



(main region), 



A\ : c > a + 5 + d 

A 2 : d > a + 6 + c, 

.D : - (a + 6 + c + d) > a, 6, c, d, 

Fi : a > 6 + c + d, 

F 2 : 6 > a + c + d. 



(1.2) 



The regions A x and A 2 are antiferromagnetic, F\ and F 2 ferromagnetic, and D 
disordered. The case d = corresponds to the six-vertex model. 

The solution of the eight- vertex model in the antiferromagnetic region is based 
on appropriate parametrization of the .R-matrix so that it satisfies some functional 
equations. We recall main points of this approach in Sec. 2. A simple local trans- 
formation makes it possible to find necessary parametrization in the ferromagnetic 
case (Sec. 3). Some nonlocal transformation ('duality') connects the ferromagnetic 
region with the disordered one (Sec. 4). A particular case of the six- vertex model 
is considered in Sec. 5. Several concluding remarks are given in Sec. 6. 

2. Antiferromagnetic Region 

In this section we recall main points of the Kyoto approach. Consider the antifer- 



romagnetic region A\. Let a\ 



(i) 



0, 1 be ground state configurations 



Vkl — v ki — (— 
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(2-1) 



Let be the set of configurations that differ from the ground state a\ l> by a 

finite number of spin flops. Only configurations from A^ U A^ contribute to the 
partition function. 
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Fig. 3 

The .R-matrix can be parametrized as follows 1 ' 7 

a(C) = a(C|«,x) = P (0 ■ ^e q 2(qx 2 )e q ,(qC 2 )e q 2(x 2 /C 2 ), 

6(0 = 6(CI«,a:) = P (0 • ^e 92 (^ 2 )e, 2 (c 2 )e 92 (^ 2 /C 2 ), 
c(C) = c(Ck,x) = p(C) ■ -e q 2(x 2 )e q ,(qc 2 )e q ,(qx 2 /c 2 ), 

1 

d(C) = d(ck, x) = /9(C) • ^e g2 (x 2 )e g2 (c 2 )e g2 (x 2 /c 2 ). 



(2.2) 



Here 



e p (z) = (*;p)oo(p/*;p)oo(p;p)oo, 

oo 



(z; P1 ,---, PN ) 00 = n (i-^r---pr)- 

ni,---,njv=0 

The function Q p (z) is connected with the standart theta-function Q(u) as follows 



0(iu) = Q q 2 (qe-™' 1 } , q = e^ 1 'I 1 ', 



with i", 7' being the standard halfperiods. 

The region A\ corresponds to the values of parameters 

0<g*<a;<C<l. (2.3) 

The limit q — > corresponds to the six-vertex model. 
The .R-matrix satisfies the Yang-Baxter equation 1 

^i?(Ci/C2) £ a 1 f 2J R(Ci/C3) e a 4 ?3^(C2/C3)^^ 

a/37 

= £ ^(C2/C3)&>(Cl/C3)?; 7 ^(Cl/C2) Q 1^. 

a'/3'7' 

(2.4) 

Grafically this equation is presented in Fig. 3. Spectral parameters Ci are attached 
to oriented lines. 
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The function p(Q in Eq. (2.2) is chosen so that the partition function 'per 
cite' is equal to unity. It is equivalent to the conditions 7 (Fig. 4) 



R(1)?H = SI'S','. 



a/3 



£1 £2 
e 4/ 



(initial condition), 
(unitarity), 



(2.5a) 
(2.5b) 



RiOlfet = R^/Oll-Z (crossing symmetry). (2.5c) 



The last condition (2.5c) is in a sense the crucial one. It is not connected with nor- 
malization. Unnormalized but obviously analytical .R-matrix p _1 (C) R(C) satisfies 
Eq. (2.5c). Therefore corresponding partition function 'per cite' p _1 (C) satisfies 
the equation 

p{x/Q = p(0- (2.5c') 

It should be mentionned that in the region Ai the crossing (2.5c) does not change 
boundary conditions (2.1). Moreover, going from ( in (2.3) to x/( we do not 
intersect the surfaces of phase transitions. Indeed, ifO<£<C<l then < x < 
x/Q < 1. But within one region the partition function per cite p _1 (C) is analytical 
and the condition (2.5c') is valid as an equation for an analytical function p(Q- It is 
just Eq. (2.5c) for the parametrization (2.2) that breaks down in the ferromagnetic 
and disordered regions. Indeed, in these regions the right hand side and the left 
hand side correspond to different phases with different greatest eigenvalues of the 
transfer matrix. The partition function per cite p _1 (C) for unnormalized analytical 
weights cannot be analytically continued over the phase transition. Hence, the 
parametrization (2.2) is useless in ferromagnetic and disordered phases. 
The solution to Eqs. (2.5) is given by 1 ' 7 



Pit) = V; <z 2 )- VC 2 ; qV(q^ 2 C 2 ; s)" 1 

(x 4 C 2 ; q, a^ooO^C" 2 ; q, x 4 ) og (qC 2 ; <?, x A ) ^{qx 2 C~ 2 ] 3, x 4 ) c 



X 



(2.6) 



(x 4 C" 2 ; q, ^ 4 )oo(x 2 C 2 ; q, x 4 )oo(?C- 2 ; <z, ^ 4 )oo(g^C 2 ; q, x 4 ) c 



We now turn to correlation functions. The Kyoto approach uses two main 
types of objects: corner transfer matrices (CTMs) and vertex operators (VOs). 

The north-west corner transfer matrix, 1 C^wiOelel •••> i s the partition function 
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Fig. 6 



on the north-west quadrant over the set with the configuration £1^2 • • • on the 
vertical boundary and e' x e'^ ■ ■ ■ on horizontal one (Fig. 5a). The definitions of the 
south-west, south-east and north-east CTMs, Cg^/(C)> Cg E (£) and C$ E ((), (Fig. 
5b) are evident. On the infinite lattice 







(2.7a) 


cf w {Q = 


-<(x/0 D "\ 


(2.7b) 


C ( s E (C) = 


:Cr ooC O-oo; 


(2.7c) 






(2.7d) 


^ne(0^se(0^sw(0^nw(0 — x2D ' 


(2.7e) 


where is an operator independent 


of C, and 










is a spin flop operator along one boundary of a CTM. Note that Eqs. 


(2.7b-d) 



follow from Eq. (2.7a) and the crossing symmetry (2.5c). It means that they only 
hold for the antiferromagnetic region. 
The vertex operator 



^ M) (0^::: = EII fl (C)JL? 

{Mfe} k=1 



fJ, =£ 



?k = -4 = (-) 



k+i+l 



k > 0. 

(2.8) 



is a partition function along a half line (Fig. 6). VOs satisfy the equations 
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e' e' 

t l t 2 



fl(Ci/Ca)&W 
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(Ci) 
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(C 2 ) = ^ 1 2 - l ' l) (C 2 )^i' 1 - l) (Ci), (2.9a) 
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1. CTMs 

2. VOs 

3. i?-matrices 



(6) 



e *?- m) (c)=*?- m) (c/*k 



(<) 
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(2.9b) 
(2.9c) 



Two first eqations follow from the Yang-Baxter equation (2.4), Eq. (2.9c) is 
evident. 

Every correlation function can be expressed in terms of vertex operator cor- 
relation functions 



F £l ... £n (Ci, • • • , Cn) = Tr (x 2D(l) *(M"0 (Ci) • • • (Cn)) , n = 0, 2, 4, • • • . 

(2.10) 

Indeed, consider, for example, the probability P ei ... £m that m spins along a row of 
parallel links take the values S\, ■ • • , e m (Fig. 7a) in the phase . Then 



p £1 - £m = Tr (cg2(c)c'iS(c)<' 1 - i 'Hc)---*g- M) (c) 

-eiei--e m \ x Ci ' ' ' ' ^C; C> ' ' ' ) 0- 



(2.11) 



= F d) 



The last line is obtained using Eqs. (2.9b) and (2.9c). In general case it is necessary 
to divide the crystal (Fig. 7b) into the CTMs 1, products of VOs 2, and a finite 
cluster of .R-matrices 3. 
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VO correlation functions can be found by solving the equations 7 

^ii-"£ n (CCl) ' ' ' ? £Cn) = F El ... e „(Cl? ' ' ' ? Cn)? 

V] R (Cj/ Cj+^SjejllF^, , ...(• • • , £j, + 1, • • •) = ^■••e,-+ie,-- G ' ' , + 1, (?» ' ' ')» 
e'e' 

■^ei£2"-e n Cl) C2? ' ' ' ? Cn) = ^62-"£m£i (C2? ' ' ' ? Cn? Cl)? 
F-l ££l ...e n (^C? C> Cl> ' ' ' ) Cn) = Fe 1 ... £n (£1, • • • , £ n ), 

e 

F w () = 1. 

(2.12) 

For example, 7 

£ £ eF^ s (c, 1) (a^r 1 ; ^)oo(xC; ^ooC-xC" 1 ; (zJooC-^-'C; 



E e i^eCC, 1) (-^C -1 ; ^ 2 )oo(-<; a^ooOcC -1 ; 9)oo(?s _1 C; 9)00 ' 

and the Baxter-Kelland formula 1 for the 'staggered' spontaneous polarization 
holds 222 2 

I \ _ ( X ? X )oo(~Qi 1)qo In io\ 

(e) -(-*V 2 )L(«;«)L" ( J 

3. Ferromagnetic Region 

Consider the eight- vertex model in the main region Ai. Consider the substitution 

m = (-) k+l »w m = (-) k+l v' kl . (3.i) 

Introduce new .R-matrix, Rp, so that 

pMfei "fc.i-l / p/ \Mfei ^fe,;-i / 7-,/ \— Mfc! ^fc,(-i q \ 

n fik-i,l"kl — \ n F) u ' v' ~ \ n F)fik-i,l-Vkl ' l^o.^aj 

^k — 1,1 kl 

or, in components, 

= c, bp = d, cp = a, dp = b. (3.2b) 

Evidently 

ap > bp, cp, dp, 

and the matrix Rp describes the model in the ferromagnetic F\ region. The 
ground state antiferromagnetic configuration a\ transforms to a ground state 
ferromagnetic configuration 



m = m = (-f (3.3) 
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under the substitution (3.1). Similarly, the set transforms to the set of 

configurations that differ from f[ by finite number of spin flops. 

The substitution (3.1) defines a natural parametrization in the ferromagnetic 
region 



a F (() 


= p(C) 


ie 92 (o: 2 )e g2 (QC 2 )e g2 (^ 2 /C 2 ), 


MO 


= p(C) 


1 

^e, 2 (o; 2 )e g2 (c 2 )e g2 (a; 2 /C 2 ), 


cf(C) 


= p(C) 


^e (?2 (^ 2 )e g2 ( g c 2 )e g2 (x 2 /c 2 ), 


d F (C) 


= p(C) 


le (?2 (^ 2 )e g2 (c 2 )e g2 (^ 2 /c 2 ), 



(3.4) 



with (/, x and £ in the region (2.3). The .R-matrix Rf(C) satisfies the Yang-Baxter 
equation (2.4) and equations (2.5a,b). The crossing symmetry is different (Fig. 8) 

R F m £ et= R F( X /O e efel- ( 3 ' 5 ) 

Then the CTMs are given by 

c { X(0 = c^ E (0 = C D(i \ 



Cf w = CW B (C) = (x/C) Di 
with some new Z)W. The VOs, <I>^(C) satisfy the equations 

£ ^(Ci/c 2 )&^^(Ci) = $S(c 2 ) *g>(Ci; 

e' e' 

^ w ^(c) = ^(c/oe DW . 



(3.6) 



(3.7) 



These equations do not relate objects with different values of i and we omit the 
superscript ^ from now on if it does not lead to a confusion. The VO correlation 
functions 

F £l ... £n (( lr --,(n)= Tr (x 2D $ £l (Ci)---$ £ „(Cn)), n = 0,1,2,-.., (3.8) 

9 



satisfy the equations 



F £l •••£„ (CCl) ' ' ' 7 £Cn) ■fei--e n (Cl> ' ' ' > Cn)) 
X] R F(Cj/Cj+l)e i eiX\F...e'.e'. +1 ---(' ' ' ,0,0+1, ' * = ^•■■£ J + i£ J •■■ (" ' * , + 1,0, * " 0, 



£'.£'-, 



-Peie 2 ---e„(^ 2 Cl! C2, ' ' ' , Cn) — F £2 ... £n£l (C2, ' ' ' , Cn, Cl), 
^ ] ^eeei---e n (^0 0, ' ' ' , Cn) = F £l ... £n (Cl> ' ' ' > Cn)> 

£ 

^() = i, F £ (i) (C) = 5 £)( _ 1)I . 

(3.9) 

The probabilities P ei ... £m are given by 

Pe 1 -s m = F £m ... £l£l ... £m (x(, C, • • • , 0- (3-10) 



The substitution (3.1) identifies Eqs. (3.9), (3.10) with Eqs. (2.12), (2.11) for even 
n. It proves the correctness of the above reasoning. In particular 

-p£i£ 2 fer £„_i£„(Cl, C2, • • • , Cn-1) Cn) = ^ie 1 e2 t - lf -e„_ie Tl (Cl! C2, ' ' ' , Cn-1, Cn)- (3-12) 

To conclude this section we shall make several remarks. 

Firstly, if one does not want to work with equations for VO correlation func- 
tions, one may use the parametrization (2.2), but with new value of the parameter 
x and new function p(Q- The connection between the parametrizations (2.2) and 
(3.4) is given in the Appendix. 

The six- vertex limit dp — > corresponds [in parametrization (3.4)!] to 

q — > 0, —q~zx = const = Xp. (3-12) 



Here xp is the parameter of the six- vertex model corresponding to x in the 
parametrization (2.2). 

The spontaneous polarization in the ferromagnetic phase is given by the same 
formula (2.13) (but with the relation (3.4) between £, q, x and weights!). In 
the six- vertex limit (3.12) we obtain (s) = 1. This result is consistent with the 
wellknown fact that the ferromagnetic phase in the six-vertex model is frozen. 

Another interesting question is the sense of VO correlation functions with odd 
n. Physically they correspond to a crystal with dislocations. It is easy to check, 
for example, that the spin at the link of a single dislocation for the eight-vertex 
model in the ferromagnetic phase is frozen (e) = 1. There is an analogue to these 
correlation functions in the antiferromagnetic region: correlation functions with 
one insertion of erL. 
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4. Disordered Region 

The duality transformation 1 

o-d = 7}( a F + b F + c F + d F ), 

i>D = \{aF + b F - c F - d F ), 

I (4-1) 
cd = ~(a F -b F + c F - d F ), 

d D = -^(of - b F - c F + d F ) 

connects the disordered region with the ferromagnetic region F\. Physically the 
duality is based on the representation of the .R-matrix 

(RlYefet = \ a »i ( R lYefet = ^l^S^, 

Weill = \b D e 1 e 3 , (lU)™ = \b D e^ (42) 

{Resist = (ReYefet = \c D E 2 E^ 

= \d D S^ (RsY s f4 = \d D E X E 2 . 

The partition function can be written as 



"k,l-l 

,l v k i 

{tfci = l,---,8} {^ki,v k i} kl 



Let us sum up over all fikh v kl- A nonzero contribution into the sum over, for 
example, fiki is given by the terms 

y {R tkl Y ki Ukj - 1 (R tk+1 y k+i ' 1 Uk+i ' i ~ i (4.3) 

Vkl 

with such tki and tk+i,i that the power of \iu is even. Now one can associate a 
dual vertex 63 64 to any R t (Fig. 9) and a dual spin to any pair (tki,tk+i 1) or 
(tki,tk,i-i) that gives a nonzero contribution: 

#1 — > ~\ ^ , R2 —> _ _, -R3 ^ _ ~\ , -R4 — >■ ^ _ , 

(4.4a) 

+ — — + — — + + 

-R5 ^ _ _|_, -R6^_|_ _ 5 -R7 — ^ _|_ _|_, -Rs ^ _ _ 
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(t,t')- 
(t, t') - 

for vertical links, and 



+ if t= 1,4,5,8, t' 
if £ = 2,3,6,7, t' 



1,4,6,7, 
2,3,5,8 



(4.4b) 



(t,f) 



+ 



if 
if 



1,3,6,8, 
2,4,5,7, 



t' 

a 



1,3,5,7, 
2,4,6,8 



(4.4c) 



for horizontal ones. Then the summation over dual spins is equivalent to the 
summation over {tki}, and this dual model corresponds to the local weights an, 
bo, cd, do- This is the Baxter's construction of duality. 1 Note that the dual lattice 
here coincides with the initial one and is not shifted by a half period as it is for 
the Ising model. 

Now we shall extend this duality to correlation functions. Let us insert into 
the lattice several spin flop operators a . If, for example, one of them is inserted 
into the vertical link kl, the expression (4.3) turns into 



Mfc-l,i v kl 



( F? ^ Mfe + 1,( ^fc + l,!-l 



fj-kl 



±Y(Rt kl Y ki Uk ' i ~ i (Rt k+1 y kA 



\Hk + l,l ^fc + lJ-1 



Mfci 



where the sign is '+' if there is no fiki in both R tkl and Rt k+1 1 and '— ' if fiki is 
there. From Eq. (4.4b) we see that the sign coincides with the value of dual spin 
at this link. It means that a spin flop operator in the initial lattice corresponds to 
a spin variable in the dual lattice. Conversely, a spin variable in the initial lattice 
corresponds to a spin flop operator in the dual lattice. It means that there is a 
relation between correlation functions in the disordered and ferromagnetic regions 

■ ■ ■ Vk m i m Vk\v x ■ ■ ■ Vk>j> n v 1 K 1 L 1 ■ ■ ■ (j1 k m l m t k' 1 l' 1 ■■■ t k' n l' n )d 

1^,! ! ! ! (4-5) 

- 2^ a kih '" (T k m l m T k' 1 l' 1 '"Tk'J'^K^ ■■■VK M L M VK' 1 L' 1 ■■■^K' N L' N ) Fl , 

i=0 



12 



where a\ { and r\ { are spin flop operators a 1 inserted into the vertical and horizon- 
tal links kl respectively. The summation over the boundary conditions i = 0, 1 is 
necessary, because the construction of duality needs summation over all configu- 
rations. 

We now turn to parametrization. The transformation (4.1) leads to the 
parametrization 

a D (0 = M P (oe q (-oe q (-x/oe q (qk)Q q (q l2 x/o/e q (-x)e q (q^x), 

MC) = Mp(0O q (0O q (x/0O q (-qk)® q (-q h x/0/O q (-x)O q (q?x), 

CD (0 = M P (oe q (-oe q (x/oe q (qk)Q q (-q h x/o/e q (x)e q (-q^x), (46) 
d D (Q = Arp(c)e g (c)e g (-x/c)e g (-^c)e 9 (^x/c)/e (? (x)e (? (-^x), 

v= 1 - 



1 (q;q)oo(x 2 ;q) 00 (qx 2 ;q) c 



X (-1; q)^ -q)^ 

The .R-matrix Rd(() satisfies Eqs. (2.4), (2.5a,b) and the same crossing symmetry 
equation as the ferromagnetic .R-matrix 

Rd(QII%=Rd(x/QI»1 (4.7) 

So the CTMs are given by 

c nw (0 = c S e(0 = C d , 

Csw(0 = C NE (0 = (x/0 D , 
and the VO, & s ((), satisfies the equations 

E^(Ci/C2)siS^(Ci)^(C2) = <MC 2 )<MCi), 



(4.8) 



-l c 2 



(4.9) 



^<M0 = <MC/£)£ D 

Therefore the VO correlation functions 

F £l ... £ „(Ci,-",Cn) = Tr (x 2D $ £l (Ci)---$ £ „(Cn)), n = 0,l,2,..., (4.10) 
satisfy the equations 

F £l ... £n (£Cl? ' ' ' ) CCn) = F £l ... £n (Cl, ■ ■ • , Cn) 
Yl R D{Cj / C,j + l)e j e j + +\F... e i, e > ] + i ...{- ■ ■ , (j, ■■■) = F... £j + l£j ...(- ■ • , (j, ■ ■ •). 

„2, 



e' e' 



F £l£2 --- £n (x Cl? C 2 ? ' ' ' j Cn) F £2 ... £n£l (C2? ' ' ' ? Cn? Cl)? 
^ ] ^eeei---e„ (^C? C? Cl? ' ' ' ? Cn) = F £l ... £n (Cl? ' ' ' ? Cn)? 



FQ = *UC) = \ 



(4.11) 
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with Eq. (3.10) for probabilities. 

Let us relate this construction to the ferromagnetic case. Let 

* e (C) = *'+(C) + (0- (4-12) 
It is easy to check using Eq. (4.1) or (4.2) that 

E^(Ci/C2)ft^yci)^(c 2 ) = ^ a (C2X(Ci), 

e' e' 

t l t 2 

Therefore $^(C) is the ferromagnetic VO. More precisely, 

^(C)~(^ 0) (C) + ^ 1) (C))fe rrom a g , (4.13) 

This imposes an additional condition on the VO correlation functions 

F ei ...e n (Cl, • • • , Cn) = F_ £l ..._ Sn (Cl, • • • , Cn)- (4.14) 

On the other hand it is straightforward to obtain Eq. (4.5) from (4.10) — (4.14). 

We ought to make a remark. We have implicitly supposed that a, 6, c, <i > 0. 
But to cover the whole region \{cld +&d +cd + <^d) > ad, &D; cd, dj} > we need 
to consider negative values of a, 6, c, gL In this case the applicability of the above 
construction can be justified by the fact that there are no phase transitions in the 
disordered region. Hence, the partition function and correlation functions may be 
analytically continued from the subregion a, 6, c, d > 0. 1 

5. Six- Vertex Model in Critical Region 

One of the most interesting problems is to obtain a descripition of the six-vertex 
model in the critical region which corresponds to the disordered region of the 
eight- vertex model with d = 0. We could set do(C) = 0, but it is not convenient, 
because the crossing transformation (4.7) does not preserve this condition. Instead, 
we apply the transformation described in Sec. 3 to the disordered region 

a' D (0=c D (C), b' D (Q=d D (0, c' D (C) = a D (0, d' D (0 = b D ((). (5.1) 

This transformation maps the disordered region onto itself. The crossing simmetry 
of the form (2.5c) 

a' D (x/Q = b' D (Q, c' D (x/C) = c' D (C), d' D (x/0 = d' D (0 

preserves the condition d' D (Q = 0. On the other hand, the equations for VO 
correlation functions and probabilities are given by (2.12), (2.11) for these weights, 
but without superscripts W, with the substitution R(Q — > R' D ((), and with 

additional condition (4.14). 
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It is convenient to use here an additive spectral parameter instead of multi- 
plicative one. Change variables 

C = x = e -nX'/2I_ (52) 

It is well known 1 that 

a : b : c : d = snh (A' — u) : snh u : snh A' : k snh A' snh u snh (A' — u'), 

where snh z = — i sn iz, and k is the module of the elliptic functions. Consider 
the limit -> -1 (k -> -1). Then 

a : b : c : d = tan(A' — u') : tan-u' : tan A' : — tan A' tanu' tan(A' — u). 

In this limit a + b — c — d = or d' D = 0. Introducing new variables 

u = 2u\ A = 2A', 

5.3) 

< U < A < 7T, 

we easily obtain 

a' D (u) = p (u) sin(A — u), 
b' D (u) = p'(u)smu, 
c D (u) = p (u) sin A, 
d' D (u) = 0, 

where the function p'(u) providing initial condition, unitarity, and crossing sym- 
metry can be adopted from Ref. 9: 

Py ] n \n) \ n) \ n J JUL r p (0)r p (A) 



p / 2p\ — u \ -p / I 2pA— u 

^ = v^^/ v ^ 

v ; r ^ (2p+l)A- M j p ^ | (2p— 1)A— 



(5.5) 



Now the equations for correlation functions in the critical region of the six-vertex 
model are evident. 

6. Conclusion 

We extended the vertex operator formalism from the antiferromagnetic region to 
the ferromagnetic and disordered ones. The crucial role belongs to the parametri- 
zation of the weights and to the crossing symmetry equation which are different 
in different regions. Baxter's symmetries prove the correctness of the result. 

In this paper we did not consider vertex operators of type II. These operators 
are very important because they diagonalize explicitly the transfer matrix and 
give the spectrum of exitations. Evidently, they can be introduced in every region 
through the connection with the elliptic algebra A qy - X (sl2)- S 
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Appendix 

In practical calculations the relations between different parametrizations may be 
useful. Here we cite them: 



Rf{C\q,x) = R{(\q,x F )\p= PF , 
Rd(C\Qi x ) = R((D\qD, x D )\ p=PD . 

Using the usual techniques of elliptic functions we obtain 

_ 

xf = —q 2 x, 

p F ((\q, x F ) = ~^-p((\q, -q*x), 
x p 

+ 7T> q D =iq, Cd = C, x D = (qzx) a , a = 



(Al) 



l°g qo logQ 2tt \ogq 

PD((p\qD,x D ) __ 1 qj ® q l( x2 D)® q l(Cp)® q l{xl/C 2 D )® q {x)e q {-q^x) 

p(C\ q ,x) ' ~Mx 2 D e q (c)e q (-x/oe q (-qk)Q q (q^x/o 

(A.2) 

Here the relation between the conjugate parameters q' and q' D is cited. 

Note that the ferromagnetic spontaneous polarization can be written in the 

form 

(qx%;qx%) 2 x (-q;q) 2 c 



(e} = 



' oo 

(-qx 2 F ; qx^lciq; q) 2 ^' 
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